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The quantization of the family of linearly polarized Gowdy T"^ spacetimes is discussed in detail, 
starting with a canonical analysis in which the true degrees of freedom are described by a scalar 
field that satisfies a Klein-Gordon type equation in a fiducial time dependent background. A time 
dependent canonical transformation, which amounts to a change of the basic (scalar) field of the 
model, brings the system to a description in terms of a Klein-Gordon equation on a background that 
is now static, although subject to a time dependent potential. The system is quantized by means of 
\^ , a natural choice of annihilation and creation operators. The quantum time evolution is considered 

. and shown to be unitary, allowing both the Schrodinger and Heisenberg pictures to be consistently 

' constructed. This has to be contrasted with previous treatments for which time evolution failed to 

^ be implementable as a unitary transformation. Possible implications for both canonical quantum 

. . gravity and quantum field theory in curved spacetime are commented. 
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I. INTRODUCTION 



> 

■ In the search for a quantum theory of gravity within the canonical approach, it has always been useful to analyze 
symmetry reduced models. On the one hand, this allows to discuss with specific examples conceptual and technical 
issues that arise when trying to conciliate gravity and quantum mechanics. On the other hand, these reduced models 

■ are usually of physical relevance in cosmology or in astrophysical situations. The most studied examples are mini- 
, superspaces , where the infinite dimensional system is reduced by symmetry considerations to a model with a finite 

number of degrees of freedom. A more interesting and far-reaching class of reduced models, where the resulting system 
is still a field theory with an infinite number of degrees of freedom like general relativity, are known as midi-superspaces 

^jTj] The simplest of all inhomogeneous midi-superspaces in pure general relativity with spatially closed spatial sections 
• • . and cosmological solutions (expanding from a big-bang singularity) is the linearly polarized Gowdy model 
. !^ This explains the considerable attention that has been paid during the last thirty years to the problem of quantizing 
^> , this model 0, 0, 1^ iSu .Ei 10- .IL] . After the first preliminary attempts to construct a quantization and obtain 
^ ' physical predictions for the Gowdy cosmologies employing conventional (but not always rigorously implemented) 
5^ ' canonical methods in quantum cosmology '3, 'B',^ , the problem was revisited using Ashtekar variables in the context 
of a nonperturbative quantization J^^. Nonetheless, it is only recently that true progress has been achieved in the 
task of introducing a consistent quantization, at least for the (sub-)model with linear polarization 9] for which the 
two spacelike Killing vector fields of the system are hypersurface orthogonal. 

The quantization proposed in Ref. 9] for the linearly polarized Gowdy model is based on the equivalence that 
exists between the set of solutions for its spacetime metric and the classical solutions for a scalar field coupled to 
gravity in 2 -I- 1 dimensions, defined in a manifold whose topology is xT^. In more detail, after a suitable (partial) 
gauge fixing of the Gowdy model, which includes the choice of an internal time, the linearly polarized Gowdy 
spacetimes are described (modulo a remaining global constraint) by a "point particle" degree of freedom and by 
a field (f) that is subject to the same equation of motion as a massless, rotationally symmetric, free scalar field that 
propagates in a fictitious two-dimensional expanding torus. The quantization of the local degrees of freedom of the 
Gowdy model can hence be confronted by constructing a quantum theory for this scalar field. The quantum Gowdy 
T'^ model is defined by introducing a representation for the field on a fiducial Fock space and imposing on it the 
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constraint that remains on the system as an operator condition, in order to finally obtain the Hilbert space of physical 
states. 

However, there is an important drawback to the quantization presented in Ref. 9]. It can be proved that the 
quantum evolution admits no implementation as a unitary transformation. Moreover, this negative result applies to 
the implementation both on the kinematical Hilbert space 10] and on the physical Hilbert space of the model JlJ. To 
make things even more subtle, it turns out that the dynamics can be approximated as close as one wants in terms of 
unitary transformations but, nonetheless, the true evolution cannot be represented by a unitary operator. Owing 
to this failure of unitarity, we do not have at our disposal a Schrodinger picture with an evolution that conserves the 
standard notion of probability [ill fl^ l . In a more pessimistic note, one might argue that such nonunitarity poses 
serious problems for a proper physical description of the system [l^l- A careful analysis of this issue would re quir e 
further discussion about the existence of the Heisenberg picture when the Schrodinger one is not available [T^. Il5j|. 
clarifying its physical validity and elucidating whether one should or not abandon the concept of unitary evolution. 
We shall not pursue this avenue here but, considering the Gowdy cosmology as a particular arena in which one is 
addressing the issue of unitarity in cosmology, we will rather show that the problems with the quantum evolution can 
be solved, at least in this case, by adopting a different quantization for the model. 

In the study of quantum cosmological models, a fundamental issue is the so-called problem of time. In general 
relativity there are no preferred foliations in spacetime and the dynamical evolution should consider all possible 
spacelike foliations. This is one of the main features of diffcomorphism invariance. Furthermore, in cases with compact 
Cauchy surfaces, dynamical evolution is pure gauge since there is no true Hamiltonian. Then, any interpretation of 
time evolution is normally obtained via a deparametrization which, in Hamiltonian language, is achieved by fixing 
the time gauge. Thus, the dynamics to be considered in these quantum cosmological models concerns the evolution 
of quantum states between Cauchy surfaces, defined by the particular choice of time gauge adopted. Different choices 
of time may lead to inequivalent quantizations. In the linearly polarized Gowdy model the system is partially 
gauge fixed at the classical level and in particular a time function t is chosen and interpreted as the time that defines 
"evolution" . The surfaces of constant t for the quantum gravity model turn out to be Cauchy surfaces of the quantum 
scalar field in a fiducial background equipped with a foliation of preferred surfaces. Furthermore, even if there is 
no preferred time in the fundamental description of quantum gravity, as well as in the cosmological Gowdy models, 
and one would only expect a genuine notion of time to arise in a certain semiclassical regime, the introduction of 
a deparametrization allows to introduce a family of true observables, the so called evolving constants of motion [T^ 
that can be associated with quantities "living at time . A clean construction and interpretation of these observables 
turn out to be possible in our case, since both a Schrodinger and a Heisenberg pictures will be shown to exist. 

This work has several aims. First, as we have just commented, we will prove that it is indeed possible to achieve a 
unitary quantum dynamics in the linearly polarized Gowdy cosmology. Therefore, no fundamental obstruction exists 
to the standard probabilistic interpretation of quantum physics in this inhomogeneous cosmological framework. An 
outline of this result was presented in Ref. |l7| . The second aim of the present paper is to systematically explore 
the canonical structure of the Gowdy model. As we will show through a detailed analysis of the implementation 
and consequences of a canonical transformation on phase space, one can arrive at a suitable field-parametrization of 
the spacetime metric of the model (i.e., to adopt an adequate choice of basic field) which allows of a fully consistent 
quantization. In addition, we also want to discuss some relevant physical phenomena that occur in the model, e.g. the 
production of "particles" by the vacuum of the cosmological system and the recovery of a time-translation invariance 
in the asymptotic region of infinite large times. 

The rest of the paper is organized as follows. Sec. reviews the standard formulation of the Gowdy model 
along the lines put forward in Refs. 0,0,0,^3 (including the realization of its nonunitary character). In Sec. IIIII 
we perform a time dependent canonical transformation and find the corresponding new Hamiltonian for the model. 
The system is then recast in Sec. IIVI as a scalar field in a static background with a potential. In that section we 
also find the general classical solution for this scalar field and the finite symplectic transformation associated with 
the time evolution. This provides the starting point for Sec. El where the quantum representation is defined and the 
Fock space is constructed. The issue of the unitarity of the evolution is analyzed in Sec. IVII It is shown that the 
symplectic transformations that define finite time evolution can indeed be implemented in a unitary way. Sec. IVIII 
studies the properties of the quantum theory, in particular its relation with previous work and the issue of "particle 
production" . Finally, we present our conclusions and some further discussion in Sec. IVIIII An appendix is added, 
where we show that the description of the model adopted in the main text corresponds to the same gauge fixing that 
had been imposed in previous treatments of the Gowdy cosmologies ■ 
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II. THE POLARIZED GOWDY MODEL 



Let us briefly review the symmetry reduced model employed in Ref. 9] to introduce a quantum theory for the 
linearly polarized Gowdy cosmologies. These cosmologies are described by vacuum spacetimes with the spatial 
topology of a three-torus that possess two commuting, axial, and hypersurface orthogonal Killing vector fields ■ They 
provide the simplest of all possible cosmological models with compact spatial sections and local degrees of freedom. 
After a gauge fixing procedure, in which all the gauge degrees of freedom are removed except for a homogenous one 
|l2L Il8| , the metric of the symmetry reduced model can be expressed in the form 

<i»' = e'e-♦'"^(-^ii^ + ^i^)»)+e-*'^iV^fa- + e*'"*<i^^ (2.1) 

n— — oo.n^O 

Here, t > is a positive time coordinate and the angular coordinates 0, a, and 5 belong to . The two Killing 
vector fields are given by da and ds, so that the metric functions are independent of these angles. The function 4> and 
its canonical momentum may depend on both t and 6. Therefore, they describe a fieldlike degree of freedom of 
the metric. All cyclic integrals in Eq. (|2.2|) are performed on the corresponding angular dependence and the prime 
stands for the derivative with respect to 6. On the other hand, Q and p are homogenous variables 112]. Moreover, p 
is a constant of motion that we impose to be strictly positive. Actually, spacetimes with p < Q can be related with 
those with p > by means of a time reversal, whereas spacetimes with vanishing p can be consistently removed (since 
the considered sector of phase space is dynamically invariant). To avoid dealing with the positivity restriction, we 
introduce the definition 

P := \np. (2.3) 

The real constant of motion P provides the momentum canonically conjugate to Q. So, the pair (Q, P) describes a 
"point particle" degree of freedom. 

To be completely rigorous, the metric (|2.1|l should include a nonvanishing 0-component of the shift vector equal 
to an arbitrary function of time Nonetheless, since this kind of shift can always be absorbed by means of a 

redefinition of the angular coordinate 9, we have obviated it. This freedom in the choice of shift appears because 
the gauge has not been totally fixed. There is still a global constraint remaining on the system, coming from the 
homogeneous part of the ^-momentum constraint: 

Co := i d9P^(j}' = 0. (2.4) 

Starting with the Einstein-Hilbert action of general relativity with AG/n = c = 1 (G and c being Newton's constant 
and the speed of light, respectively) and after the above reduction process, one arrives at the following action for the 
model (modulo the constraint Cq — and spurious surface terms): 



Sr^ dt[PQ+ (tde 



ti 



P4><P - 



nr = ^^[Pl + t\cj^'f]. (2.5) 



where the dot denotes the derivative with respect to t. The reduced Hamiltonian is 

Hr = <j> dOT-Lr. (2.6) 

This Hamiltonian does not depend on the "point particle" degrees of freedom. Hence, Q and P are constants of 
motion and a nontrivial evolution occurs only in the field sector of the system. In addition note that, were not for the 
explicit time dependence of Ti^, the reduced Hamiltonian would be that of a massless scalar field with axial symmetry 
in a static background. 

The (reduced) phase space of the system can be decomposed as a direct sum, = Fq © F, where Fq and F 
contain the "point particle" and the fieldlike degrees of freedom, respectively. They admit as coordinates the canonical 
pairs {Q,P) and (0, P^). Owing to the presence of the global constraint Cq = 0, the space of physical states does 
not really correspond to F^, but rather to a submanifold of it. However, since this submanifold is nonlinear, the 
reduction by the constraint is postponed to the quantum theory, where it is imposed as an operator condition on the 
(kinematical) quantum states. 
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The Hamiltonian equations derived from Hr require that the field (j) and its momentum satisfy 

^=^, P^^t4>". (2.7) 
Combining them, one concludes that is subject to a wave equation: 

0+^-</>" = O. (2.8) 
We will call tp any smooth solution to this equation. Let us define fn{t^ :— fn{t) exp[in9] for all n G Z with 

fo[t) ■= 7=^, Jn(t) ■■= 1=— if ntO. (2.9 

V47r V8 



Here, Hq is the zeroth-order Hankel function of the second kind 1^]. We can then express all solutions Lp in the 
generic form 

oo 

p{t,9)^ J2 [A,Jn{t,9) + A:j:{t,e)]. (2.10) 

n— — oo 

The symbol * represents complex conjugation, and the A„'s are complex constant coefficients. These coefficients must 
decrease faster than the inverse of any polynomial in n as |n| — ^ oo in order to guarantee the pointwise convergence 
of the series (|2.10|) . 

From action 1)2.5(1. one obtains the following symplectic structure on the space {p} of smooth solutions to Eq. 



Q{ipi,ip2) ^ j) d9 [ip2tdtipi - ipitdtip2] , (2.11) 

where the integral is taken over any t =constant slice. The set of mode solutions {fn{t, 9), fn{t, 9)} (with n G Z) is 
complete and "orthonormal" in the product (//, /„)^ = -zri(/;*, /„), in the sense that {fufn)ip = Sin, {fi,fn)v = -<5/n, 
and {fi,fn)ip — 0- As a consequence, it is not difficult to check that the complex conjugate constants {An, A* } behave 
as pairs of annihilationlike and creationlike variables under this symplectic structure. 

Returning to Eq. (|2.8|l . it is worth pointing out that it is formally identical to the Klein-Gordon equation of a 
free massless scalar field propagating in a fictitious three-dimensional background {A4 ~ x T^,g^^^), where the 

background metric is g'^^ — —dtadti, + d9ad9i, + t^daadab, again with t g R+ and 9,a G S^. We can then identify f 
(the field part of phase space) with the canonical phase space of the scalar field in that background, while the space of 
smooth solutions can be considered as the covariant phase space of such a Klein-Gordon field. That is, cj) and can 
be viewed as the configuration and momentum on the constant-time section ~ of the scalar field (p propagating 
in {M,g^^'>). 

A consequence of this equivalence between the gauge fixed Gowdy model and a Klein-Gordon field is that, in this 
description, the problem of quantization of the local degrees of freedom reduces to the construction of a quantum 
theory for the axially symmetric massless scalar field p in the background (Al,g^^''), which is an expanding torus. 
Employing this fact, the quantum Gowdy model was defined in Ref. Q by introducing a Fock representation for 
tp and, on the fiducial Fock space obtained in this way, imposing the global constraint (|2.4|) as an operator condition 
in order to get the physical Hilbert space. More precisely, taking into account the field decomposition (|2.10l) [and 
remembering definition l|2.9|) ]. the symplectic vector space S :~ (fi, {(^}) can be endowed with the f2-compatible 
complex structure J : S ^ S: 

J[fnit)]=lfnit), J[m]=-^fm. (2.12) 

Using this complex structure, it is straightforward to construct from 5* the "one-particle" Hilbert space of the theory, 
Ti.. This space allows to define in turn the (symmetric) Fock space !F{'H) on which one can introduce the formal 
field operator, expressed in terms of annihilation and creation operators that correspond to the positive and negative 
frequency parts determined by the complex structure J. Finally, one specifies the explicit operator that represents 
the constraint Co quantum mechanically on !F{T-L). The physical Hilbert space .?>hys is supplied by the kernel of this 
operator. 

In the quantization sketched above, however, it is known that the dynamics dictated by the Hamiltonian Hr 
[see Eqs. (|2.5f) and H2.6|l ] cannot be implemented as a unitary transformation. This result applies not only to the 
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kinematical Fock space ^Fifi.) |3, but also to the physical Hilbert space ^phys The antilinear part of the 

Bogoliubov transformation that implements the classical dynamics in the quantum theory, providing the relation 
between annihilation and creation operators at different times Iq and t/, has a single contribution for each of the 
field modes. The corresponding Bogoliubov coefficient can be deduced employing the "orthonormalization" of the 
set \fn (t, 0), f^(t, 9)} (with n e Z) and formula p.llfl . For the nonzero modes {n ^ 0), the coefficient is given by 

D^{tf,to) = i2n[f*Mtfdtf:Atf)~mtf)todtf:{to)] (2.13) 
= ^[iI*(|nK/)toi?i*(|nKo)-i/o(l"|io)i/i?r(l"|i/)], (2.14) 

where Hi is the first-order Hankel function of the second kind Since the sequenc e ipn( tf,t())} is not square 

summable for generic positive times tf to (IQj , unitary implementability is impossible |2Ct l2ll | . 

Actually, the failure of square summability can be easily derived by making use of Hankel's asymptotic expansions 
for the functions Ho and Hi Up to corrections of relative order 1/t in these expansions. 



V '^fI* V 



From these formulas and Eq. (|2.14|) . one obtains at leading order 



|2 _ (^/ ~ ^o)^ 



|A.(t/,to)r«^^^^. (2.16) 

Thus, for asymptotically large values of \n\, |D„(t^,io)P differs from zero if tf and to do not coincide, so that the 
sequence {-D„(t/,to)} is not square summable. 



III. NEW FIELD-PARAMETRIZATION 



The above discussion shows that the failure of unitary implementability of the dynamics can be blamed on the 
inadequate asymptotic behavior of the Bogoliubov coefficient Dn{tf,to)- Furthermore, a close look at Eq. (|2.13|l 
reveals that this problematic behavior can be traced back to the appearance of the factor t in the symplectic structure 
(12.11(1 . The obvious way to get rid of this unwanted factor in the symplectic (two-)form is to absorb its square root, 
^/t, in (p. Notice that such rescaled solutions are nothing but the classical smooth solutions for the rescaled field 
The proposed change results then in a rescaling of the complete set of mode solutions {fn{t, 0), f*{t, 6)} [n G Z) 
into the set {g,i(t, 9),g^{t, 6)} := {Vtfn{t, 6), Vtfn{t, ^)}i which is again a complete set of solutions, but now for the 
classical equation of motion satisfied by the rescaled field. 

Actually, the asymptotic behavior of the nonzero mode solutions fn{t,6) strongly suggests the proposed rescaling 
also from the following related perspective. Using Hankel's expansions for the nonzero modes we see that, both for 
asymptotically large wave numbers \n\ and for asymptotically large times t, 

pin/4 

/„(i,0)«^==e-(l"l*-ne). (3.1) 

Therefore, in these asymptotic regimes (and modulo negligible corrections compared with the unity) the functions 
Vifnit, 0) (with n ^ 0) behave like the standard nonzero mode solutions corresponding to an axially symmetric, free, 
and massless scalar field propagating in the static background (tW, g'-'^-'), with M ~ x and 

5afc = -dtadtb + d9ad9b + daadab- (3.2) 

Thus, the change from Lp to \Jtif> should provide us with a description that corresponds asymptotically to a massless 
scalar field in Minkowski spacetime (except for the topology). We recall that time evolution between any two fiat 
Cauchy surfaces is unitarily implementable for the latter system. 

With these motivations, we will now proceed to reformulate the linearly polarized Gowdy model by considering as 
our new covariant phase space the set {\/t<~p]. We will also see that, for asymptotically large values of the dynamics 
of the (nonzero modes of the) system is indeed dictated by the Hamiltonian of a collection of harmonic oscillators 
with frequencies a;„ = l^^l. 
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Let us hence start by multiplying the field (/) by the factor \fi and completing this rescaling into a time dependent 
canonical transformation in the most straightforward way^ , namely 

X := Px •= ^Pl>- (3-3) 

Obviously, the canonical variables (Q, for the "point particle" degrees of freedom need not be changed. The above 
transformation allows us to rewrite the Hamiltonian H2.6() as that of an axially symmetric free field, 

Hr = l<fd9[P^ + ix'r]- (3.4) 
On the other hand, given the periodicity of the system in 9^ we can expand x Px Fourier series, 



X 

%n6 



oo 



n— — OO ^ n— — oo 



It is not difficult to check that these (imphcitly time dependent) Fourier coefficients form canonical pairs, with 
being the momentum conjugate to X(n)- In terms of them, the Hamiltonian Hr adopts the expression 

^ CO 

^'■=2 51 [4"^^i""^+^'x(n)X(-„)]- (3-6) 



n— — C30 



Then, Hr can be equivalently interpreted as describing a free particle (the zero mode n = 0) and a combination of 
harmonic oscillators with frequencies equal to |n| [two for each value of n 7^ 0, corresponding to the real and imaginary 

parts of X{n) ^-nd "']. Our mode decomposition leads in this way to a natural choice of annihilationlike variables 
(up to trivial linear combinations), 

\n\X(n)+iPx^ ,„ „s 

On = (3.7) 



for all n ^ 0, with creationlike variables obtained by complex conjugation. 

However, since the change ()3.3|l is a time dependent canonical transformation, the dynamical evolution in the 
new field-parametrization is not generated by Hr anymore. Instead, the new reduced Hamiltonian of the system 
is Hr = Hr + dtF , where the partial derivative refers only to the explicit time dependence and is a generating 
functional for the canonical transformation. For instance, we can choose F = - / F^x/V^- Then dtP ^ § P^x/C^t) 
and 



2 I /,/\2 



PxX 



(3.8) 



One may decompose Hr in Fourier modes in a similar way to what we did with Hr and express the result in terms of 
the annihilationlike and creationlike variables introduced in Eq. 1)3.7(1 . Modulo a point particle, the new Hamiltonian 
corresponds to an infinite number of standard harmonic oscillators in the limit i — > 00, as expected from Eq. (|3.8|) . 

In spite of this good feature, there is a serious problem that makes us disregard Hr as a suitable Hamiltonian for 
the Gowdy model. In the Fock representation that (for the nonzero modes of the system) determines the choice H3.7|) 
of annihilationlike and creationlike variables, one can easily check that the vacuum does not belong to the domain 
of the operator counterpart of Hr'. the action of the operator on the vacuum has infinite norm. As a consequence, 
it follows that such an operator cannot be defined on the dense subspace of the (kinematical) Hilbert space formed 
by the states with a finite number of particles (indeed, none of these states has a normalizable image) . Actually, this 
problem appears owing to the presence of the cross term § P^xl i'^t) in the Hamiltonian. Nonetheless, this term is 
negligible in the limit of infinite large times and is not required to arrive at the desired asymptotic behavior for the 
system. 



^ This kind of transformation was already considered in Ref. |3 although in a different and restricted context, namely the study of the 
WKB regime. 



7 



In fact, the commented cross term can be eliminated by including a linear contribution of the field x in its canonical 
momentum. As a result, the time dependent canonical transformation (|3.3|) is replaced with the following one: 



1 



This transformation is generated by the functional 

F ^- j de 

Then, the reduced action and Hamiltonian become 

Sr = J^' dt (^PQ + 



Vt it 



d9 



dO Hr 



-(tde 



4^2 



(3.9) 

(3.10) 

(3.11) 
(3.12) 



We notice that the new reduced Hamiltonian is just that of an axially symmetric Klein-Gordon field propagating in 
fact in the fictitious static background {A4 ~ x T'^^g^^^) [see Eq. (|3.2|l ]. though now subject to a time dependent 
potential that corresponds to an effective mass equal to l/{2t). Note nevertheless that this potential vanishes in the 
limit of large times. 

It is worth remarking that the time dependent canonical transformation H3.9|) . introduced to recast the symmetry 
reduced model in terms of a new canonical set of variables (CiPj); amounts just to a field-reparametrization of the 
spacetime metric of the linearly polarized Gowdy model. The construction of the reduced model in terms of the 
canonical pair (5, Pj) is completely parallel to that explained in Ref. jl^ for the variables (0, P^) (which is essentially 
the description considered in Refs. |5ll^ lTolflH ). the only difference being the distinct parametrization of the metric. 
In all other respects, the gauge fixing and reduction process is the same, including the conditions imposed to fix 
almost entirely the gauge freedom. We show this in detail in Appendix IXI 



IV. DYNAMICS IN THE NEW DESCRIPTION 

Varying action H3.11|) one recovers the result that the "point particle" degrees of freedom (Q, P) remain constant 
in the evolution, whereas the dynamics in the field sector is dictated by 



The field ^ must then satisfy the second-order differential equation 

From Eqs. (|2.10|1 and (j3.9|l we have that all smooth solutions Q to Eq. (j4.2|) have the general form 

oo 

at,e)^ J2 [Angn{t,e) + A*^g:{t,e)]. 



(4.1) 



(4.2) 



(4.3) 



We recall that gn{t,9) := ^/tfn{t,6). As anticipated at the beginning of Sec. Ill, the set of mode solutions 
{gnit,9),g*^{t,d)} is complete. 

On the other hand, the symplectic structure on the field sector of the canonical phase space is 



A(Ki, PcJ, [6, P«J) = fd9 (6P6 - CiPc J . 



(4.4) 



Taking into account the first of Eqs. (|4.1|l . this leads to the following symplectic structure on the space {(} of smooth 
solutions: 

^{(1X2) = {d9{<:2dtCi - CidtG) ■ (4.5) 
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Comparing this with the symplectic structure (|2.11(l for the "old" field-parametrization of the Gowdy model, we see 
that the problematic factor t has indeed disappeared. In addition, the set of mode solutions {gn{t,6), g'^{t,9)} is 
"orthonormal" in the product {gi,gn)c = ~i^{9i ^9n), as required for consistency since the set {fn{t,0), fn{t,d)} is 
"orthonormal" with respect to the corresponding inner product in {fl,{ip}). 

Therefore, for the reformulated Gowdy model, the field sector of the covariant phase space is the symplectic vector 
space S := (51, {C}), which can be coordinatized either by the set {(,} of smooth solutions to Eq. 14. 2|) or by the 
complex constants of motion {Ao, {Am}), where Aq :— {Aq, Aq) and Am '■= {AmTA'^^A-m, ^-m) foi' all m e N — {0}. 
Remember that these constants are of the annihilation and creation type. Alternatively, one can consider the canonical 
phase space of the model, whose field sector is the symplectic vector space F :— (A, {(^, P^)}), coordinatized by the 
configuration ^ and momentum of the massive scalar field. Expanding in Fourier series our canonical variables ^ 
and P^, as we did for x and in Eq. H3.5() . we can equivalently adopt as coordinates the set of (complex) canonical 

pairs P^ ""*}, with n G Z. We emphasize that these Fourier coefficients depend implicitly on the time coordinate 

t. 

Analogously to our definition in Eq. H3.7|l . we now introduce annihilationlike and creationlike variables for the 
nonzero modes n 0, 



For convenience, we introduce a similar change of variables for the zero mode, although the behavior of this mode 
corresponds to a free particle in the limit of asymptotically large times, rather than to an oscillator: 

bo = ^ ^ , K = ' ^ ^ ^ ■ (4.7) 

These transformations are canonical inasmuch as the pairs 6„ and i6* are canonically conjugate for all values of n. 
Hence, we have that F can be alternatively coordinatized by the complex variables {Bq, {Bm}), where Bq := (6oj^o) 
and Bm ■= {bm, b*^,b-m, b*_m) for all m G N — {0}. The physical phase space consists then of those states (Bo, {Bm}) 
in F that satisfy the constraint Co, which can be expressed in the form 

C30 

Co=J2m {b^bm - b*_„,b-m) = 0. (4.8) 

m— 1 

Since this constraint defines a non- linear submanifold of F, we leave the corresponding reduction to the quantum 
theory. 

The dynamics on F is dictated by the reduced Hamiltonian l|3.12|l . which in terms of the new set of variables is 
given by 

Hr = ^ [uj(m,t){bmbm + b*_,J}-m) + P(m.t)(C&-m + bmb-m)] , (4.9) 
m=0 

where 

-(o.*)-i + I^' Pio,ty.= -l + j^, (4.10) 

for the zero and nonzero modes, respectively. In coordinates (C(o)7 -Fj"'')! the zero mode part of this Hamiltonian //r"'' 
can be written 

In agreement with our above comments, this Hamiltonian describes a harmonic oscillator with unit mass and frequency 
uj = l/(2<:) that behaves asymptotically like a free particle. On the other hand, for the nonzero modes, we get 
limi_,oo '^(m,t) = and limt_,oo P(m.t) = C Thus, asymptotically, the dynamics of the nonzero modes corresponds in 
fact to that of a collection of harmonic oscillators with frequency ojn = \n\. Besides, in contrast to the situation found 
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in Sec. Ill with the field x, the vacuum will be in the domain of the Hamiltonian operator since P(m,t) turns out to 
be square summable. 

The map from the covariant phase space S to the canonical phase space F is given in the case of the zero mode by 
Bo(i) = Wo(t)Ao (treating Bq and Ao as row vectors), where 



Woit) 



roit) so{t) 



with 



2t 



roit) := V^goit,e) 1 + + so{t) := V^glit.O) U + - 



2t 



(4.13) 



(4.14) 



For the remaining modes (m e N — {0}) the map is Bm{t) = W{xm)Arm where we have defined x„ 



( c{Xrn) 




V d*{Xm) 



and 



c(^m) 

d{x„i) := 



d{x^) 

C*{Xra) d*{Xm) 

d{Xm) c{x,n) 

C*{Xm) 



HoiXm) - iHi{Xm) 
HoiXrn) - iH1{Xm) 



(4.15) 



i 



Since 



ko(i)l' - ko(OI' = 1, Hx^)? - \d{xm)? = 1, 



(4.16) 
(4.17) 

(4.18) 



for all t > and m € N — {0}, the maps Wa{t) and W{xm) are Bogoliubov transformations. Hence, the map from 
S" to F is a time dependent canonical transformation. A generating functional for this transformation (that depends 
on some appropriately chosen complete sets of compatible components -under Poisson brackets- both for S and F) is 
F = EmsN ^™ with 



Fm{t) = [c{Xjn)A^yn + d{Xra)A. 



- ib,n [d [xm)A—m + c (a;„j)j4^ 



TO 7^ 0. 



(4.19) 
(4.20) 



A straightforward calculation shows that dtF — Hr- Since the evolution in S is frozen, an initial state 
(Bo(io): {l3m{to)}) in F at time will evolve to a state (Bo(i), {Bm{t)}) at time t according to 



Bo(t) - Wo{t)Wo{to)-^Bo{to), Brn{t) = W{x^)W{xl^)-^B^{to) 



(4.21) 



where x^ :— mto. In other words, the transformation H4.21(l is the integral curve of the Hamiltonian vector field 
— A^^VsHr on F, with end points at {Bo{to), {Bm{to)}) and {Bo{t),{B,n{t)}). Alternatively, it can also be 
viewed as the map that relates copies of F at different times, e.g. {{Bo{to),{Bm{to)})} at to with {{Bo{t),{Bm{t)})} 
at t. 



V. QUANTUM TIME EVOLUTION 

Given a Cauchy surface in (Ai ~ M+ x T^, g'-'^''), for instance the surface t — to, one obtains a one-to-one correspon- 
dence between the spaces S and F by means of the maps Wo{to) and W{mto) [see Eqs. ()4.13|l and H4.15|l ]. Employing 
the inverse of these maps at t = to, one can then rewrite expression l|4.3(l in terms of a new set of "orthonormal" mode 
solutions {Gn{t, 0),Gl{t, 9)), 

OO 

C{t,e)^ [Gn{t,0)bM + GUt,0)bn{to)]- (5.1) 



n— — OO 
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For the zero and nonzero (n 7^ 0) modes, respectively, these new mode solutions are given by 

Go(i,0) = VtK{to)Mt,9) - s*{to)f*{t,9)], (5.2) 
Gnit,e) = )i?o(x|„|)-d*K|)i?o*(^|«|)le"^ (5.3) 



Here, x\n\ = \n\t and x°^^ = \n\to. 

Instead of t = to, we could have considered a different Cauchy surface t = T. In such a case, the field ({t, 9) would 
have adopted an expression similar to (|5.1() . but now in terms of the set of coeflicients {6„(r), 6* (T)} and the "orthonor- 
mal" mode solutions, {G''n\t^ 6), ci"^^* {t, 9)} that are obtained by replacing to with T in Eqs. (|5.2I) - (|5.3|I . Therefore, 
associated with a uniparametric family (UF) of Cauchy surfaces T G [to,tf] C M+, there exists a UF of "orthonormal" 

mode solutions {G^n"\t,9),G^n^^* {t,9)}Te[to,tf]i as well as a UF of copies of F, namely {(Bo(T), {^^(T)})}^^^^^ 
which are related via the evolution map (|4.21l) . 

Denoting Gn(t) :— G„(t, 0) exp[— m0], the explicit decomposition of the solutions in complex conjugate pairs pro- 
vided by Eq. H5.1|) allows one to introduce the following Jl-compatible complex structure: 

J[G„(<)] - *G„(0, J[G:(<)] = -tGUt). (5.4) 

Given a UF of Cauchy surfaces, we will also get a UF of fi-compatible complex structures Jt, namely those defined 
by the UF of "orthonormal" mode solutions {Gn {t,9),Gn (i, ^)}T6[to.t/]- 

MGiPit)] = ^&P{tl JT[G^P*{t)\ = -zGf )*(t), (5.5) 

where again G^\t) Gli^-'(t, 6*) exp[— m6']. Thus, for each copy |(Bo(r), {,Bm(r)} )} of F, we obtain a natural 
complex structure Jt '■ S ^ S . Since the copies of F are related by the evolution map H4.21|l . the UF {JT}T£[to.tj\ 

just the set of complex structures induced by time evolution. Obviously, J = JT\T=to and Gn{t,9) — Gn {t,9)\T=to- 
Starting with {S,J), we can construct the "one-particle" Hilbert space TC. It is the (Cauchy) completion of the 
space of "positive frequency" solutions 5"+ := {C^ = {(^ — iJ()/2} with respect to the norm ||C^|| = \/ (C^X^)- Here, 
(•,•) denotes the Klein-Gordon inner product: {C~^,^^) = — ifl(C~,T'^) where (~ G H (i.e., the complex conjugate 
space of Ti.). The (kinematical) Hilbert space of the quantum theory is the symmetric Fock space T{Ti) constructed 
from the "one-particle" Hilbert space. That is, 

HU) = ®T=o , (5.6) 

where ^^^^Ti. is the Hilbert space of all fcth rank symmetric tensors over H.. Following this prescription, we can write 

the formal field operator C, in terms of annihilation and creation operators corresponding to the positive and negative 
frequency decomposition defined by the complex structure J: 



00 

C(t;0)- [Gn{t,9)bn + Gl{t,9p^. (5.7) 



n— — C30 



Remembering expression (I5.1|) for the classical solutions, we see that we might have obtained this field operator by a 
straightforward assignation of operators to constants of motion. This is the Schrodinger picture, where the complex 
constant coefficients {&ri(to), feJi(to)} are promoted to annihilation and creation operators {hn{to) = bn,bn{to) = bl^}. 
We note that {{S, JT)}Tm+ leads to the UF of Fock representations |(J^(7^t), {b„{T), bl{T)})\ . Clearly, the 

Fock representation constructed from (S, J) belongs to this family and corresponds to T = tg- 

On the other hand, in the Heisenberg picture, time evolution for operators is determined by the Bogoliubov 
transformation H4.21|l . Introducing then the notation bi^^to) :— bn, one obtains the following relation between 
the annihilation and creation operators at Iq and a different time t: 

bi"\t) = an{t,to)bi"\to) + Pn{t,to)U"J\to), (5.8) 
where the Bogoliubov coefficients for the zero modes are [see Eq. (|4.14l) ] 

ao{t,to) = ro(t)ro(to) - so{t)so{to), f3o{t,to) = so{t)ro{to) - ro{t)so{to), (5.9) 
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while for the rest of modes with n G Z — {0} one gets [see Eqs. (|4.16ll and (I4.17() ] 

Unit, to) = c(a:|„|)c*(a;|'„|) -d(x|„|)d*(a;|^„|), /3„(t,io) = rf(a;|„|)c(x|'„|) - c(x|„|)rf(a:|^„|). (5.10) 

Interchanging the roles of t and to in Eq. I|5.8|l . we can write the annihilation operator at time to in terms of 
the annihilation and creation operators at time t. By substituting the result in Eq. H5.8I) . and using the relations 
|ap — |/3p = 1, a„(t,to) = a*(to,t), and /3„(t, to) — —Pn{ta,t), we then get that the creation operator at time to is 
equal to l3^{to,t)VL^ {t) + a;(to, t)6;r^' (t). Hence, the field operator can be fully expressed in terms of operators in 
the Heisenberg picture: 

C(t; 9)^-^ (r^(t) + hi"^\t)) + Y: -7=Q {^'"-'bi^Ht) + e-%(^^Ht)) ■ (5.11) 

VI. UNITARITY OF THE EVOLUTION 

The time evolution in the Heisenberg picture described in the previous section is unitarily implementable on 
the (kincmatical) Fock space T(T-L ) co nstructed from {S,J) if and only if the sequence {/3„(t, to)} that appears in 
relation (|5.8() is square summable |2C| (see also Ref. |2l|). Let us remark that unitary implementability amounts 
to unitary equivalence between all the Fock representations in the UF under consideration. Since two complex 
structures Jti and lead to unitary equivalent representations of the canonical commutation relations if and only 
if their difference (Jtj — Jt^) defines a Hilbert- Schmidt (HS) operator, either on or TIt^ (see e.g. 0|), we have 
unitary implementability if and only if Jt := (J — Jt) is HS for every T € K"*". In fact, it is not difficult to see that 
Jt is HS if and only if the sequence {/3„(T, to)} is square summable (as it should be, because the requirement of 
unitary equivalence between Fock representations is just a reformulation of the unitary implementability condition). 

Let us discuss then the square summability of the sequence {/3„(t, to)}. Since /3n(t, to) = /3_„(t, to) and, in addition, 
summability does not depend on the contribution of a single term (e.g. n = 0), it suffices to analyze the sequence 
{/3,„(t,to)} with m e N — {0}. We start by showing that the sequence {(i(mt)} (m e N — {0}) is square summable for 
all t > 0. From the asymptotic expansions of the Hankel functions for large positive arguments, we know that jl9| | 

Ho{x) = \[^[P{Q.x)~iQ{Q,x)]e-'-e^^'\ (6.1) 

V -KX 

Hi{x) - \f^[Q{l,x)+tP{l,x)]e-'^e'^/\ (6.2) 



where 



P(^.,x) -.^ l+p(^,x)^l + f;(-l)'=|^, (6.3) 



k=l 



Here, v = Q^l and {y, m) are Hankel's symbols, which in terms of the gamma function are 

Substituting Hankel's expansions in Eq. H4.17|l . the square modulus of d[xm) in the asymptotic region Xm 3> 1 
becomes 

p{{),Xm) - p(l,Xm) Q(0,Xm)\^ , ( Q{0, Xm) - Q(.l, Xm) , l+p(0,Xm)^'^ 



=[ 2 i^j +1 2 + ixm ) ■ ^'-'^ 

Employing the expressions of the functions p{v,x) and Q{v,x), it is a simple exercise to check that the last term 
in roimd brackets is o{l/x%^) at infinity^, whereas the first one presents the behavior l/{Axm)'^ + o(l/x^). So, for 



^ We say that a function f(x) is o{l/a;") when x — > oo if the product x"f{x) tends to zero in this limit. 
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» 1 one gets 



This implies, for instance, that one can find a positive constant C so that < l/{ix„iY if > C [because 

1/(3)"* > 1/(4)'']. For every given value of t > 0, let us call Mq := int(C/i) < oo, where int(x) denotes the integer 
part of X. We then obtain 

oo Mo -J oo 

m=l m=l ^ ' m=Mo + l 

In the last inequality we have used that the first sum involves only a finite number of (well defined and bounded) 
terms and that the sequence {1/m^} is summablc. This proves that the sequence {d(mt)} with m G N — {0} is square 
summable for all positive times t > Q. 

Given two values of the time coordinate, t > Q and to > 0, we know from the square summability of the sequences 
{d{mty\ and {(i(mto)} [and without appealing to the explicit form of d(a;)] that there exist two integers mo and mo 
such that |(i(mt)| < 1 for all m > mo and \d{mtQ)\ < 1 for all m > mo- Remembering that |cp = 1 + jdp, one also 
has that |c(mt)| and |c(mto)| are smaller than \/2 for all m > Mi := max(mo,mo). In this case, we obtain from Eq. 



Using the inequality 



we then conclude that 



\Pmit,to)\ < \/2 (|d(x„)| + M(xOjl) . (6.9) 
{\dix,^)\ + \dixl)\f < 2 {\dix^r + , (6.10) 



Ml 



E IAn(t,MP< E l/5™(^'io)P+4 E + M(^?„)P) • (6.11) 

m—1 m—1 m—AIi + l 

Provided that |/3m(t,to)| is bounded from above for all t,to > and m e N — {0}, the first sum in the right hand 
side of Eq. (|6.11|) is bounded since it contains only a finite number of contributions. The second sum is bounded 
from above as well, because both sequences {d{xm)} and {d{x^)} are square summable. Therefore, the sequence 
{l3rn(t,to)} is square summable for all t,to > 0. 

As an aside, let us comment that Hankel's expansions H6.1|l and H6.2|l provide the following asymptotic behavior of 
the Bogoliubov coefficient /3m(t,to) for large m: 



Pm{t, to) — — J. 

16m^ 



cos(mAt) — i ( ^ + ) sin(mAi) 



{t-Atf J V(i-A<)2 i2 



o[^], (6.12) 



m^ 



where At := t — t^. Obviously, in agreement with our above comments and in contrast to the situation found in Eq. 
()2.16(l . the dominant term in the asymptotic regime is square summable. Moreover, it is not difficult to sec that, 
when t >> 1 and At/t << 1, the above expression continues to be valid for all m > if one merely replaces o(l/m2) 
with 0(1/^2). Approximating to the same order the rest of terms in the expression, we arrive at a square modulus for 
Pmit, to) with the asymptotic form 

1^ , ,,9 sin^(mAt) /l\ , s 

\Pm{t,to)\'^ ' +o[-j]. (6.13) 



64(mt)4 

For asymptotic large times we see that, if we keep At fixed, the square modulus of /3m (t, to) decreases to zero as 1/i*. 
Since the number of "particles" produced by the vacuum in the nonzero modes is given by the sum of the sequence 
{\Pm{t, ^0)1^)7 in a fixed lapse of time At the "particle" production will be attenuated as time increases. In particular, 
regardless of the fixed value of At, one gets the bound 

E ^o)P < E + o (^) . J^Z(4) + o (^) , (6.14) 

m—1 m—1 ^ ^ ^ ^ 
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where Z{4:) = 7r^/90 is the Riemann function Z{j) at j = 4. We will discuss the issue of "particle" production in 
more detail in the next section. 

The proved square summability of the sequence {Pmit,to)} for all positive times t and to ensures that the time 
evolution is unitarily implementable on the (kinematical) Fock space !F{'H) [and that the Fock representations of 
the introduced UF are all unitarily equivalent], so that probability is preserved. Moreover, one can check that the 
evolution H5.8|l leaves invariant the constraint 

oo 

Co ^ ^ TO {blJjjn - V_,J)-ry?j , (6.15) 
m— 1 

which imposes the condition that the total momentum of the field ^ in the 0-direction be equal to zero [see Eq. H4.8|) ]. 
This invariance guarantees that the dynamics is unitarily implementable not just on TiTi), but furthermore on the 
Hilbert space J-phys(7i) of physical states, which are the states that belong to the kernel of the constraint Cq. In 
conclusion, we have shown that the quantization put forward for the polarized Gowdy model is such that the 
physical evolution is unitary. 



VII. FEATURES AND CONSEQUENCES OF THE QUANTUM EVOLUTION 

In this section, we want to clarify certain mathematical aspects of the quantization and evolution proposed for the 
Gowdy model and discuss some of their physical consequences, including the cosmological production of "particles" 
by the vacuum of the theory. We divide this analysis in several parts. 

(1) Associated with the field decomposition (|4.3|) we have the il-compatible complex structure 

Jj5„(t)] =i5„(i), Ja[9*n{t)\^-i9n{t), (7-1) 

where 5„(t) '■— gn{t, 9) exp[— m6']. Starting with (il, {C}, Jg) but adopting (Ag, {-4,„}) as coordinates for {C} instead of 
(Bo(io), {'Bm(^o)}), one can construct the "one-particle" Hilbert space Tig as well as the corresponding symmetric Fock 
space TiTig). This Fock space would now provide the (kinematical) Hilbert space of the quantum theory. Defining 
J := J — Jg,we have that 

\ £ {G^,J^J[G^a])H = |.o(io)|' + 2f; |d(x^)|2 <oo. (7.2) 

n— — oo m— 1 

We can therefore assure that {T{l-Lg),{Am A]^}) and (JF(7^), {6„, foj^}) are unitarily equivalent Fock representations. 
Besides, since the sequence {d{xm.)} is square summable for all i > 0, the unitary equivalence holds regardless of the 
value chosen for the instant io > in the construction of the representation (J^('H),{6„ = bn{to),bn = ^n(*o)})- In 
this sense, the role of the time of reference is irrelevant. 

In the considered description, on the other hand, an initial state (An, {^m}) at time to evolves to the final state 
{Ao,{Am}) at time t according to 

Ao = Wa\to)Wo{t)Ao, Am = W-\x"jW{x^)A„,. (7.3) 

The antilinear part of the map 1)7. 3|l is obviously square summable and, consequently, the classical dynamics is 
unitarily implementable with respect to the Fock representation {TiTig), AJ^}), as required for consistency with 
the unitary equivalence between this representation and {J-{Ti,), 6jj}). 

However, it is worth pointing out that the transformation (|7.3|l does not represent the total change in time, which 
is actually dictated by K :— Hr — dtF = [see Eqs. (|4.19|l and (|4.20|l ]. It is rather the relation between constants of 
motion that generates the Hamiltonian 1)4. 9|l . written in coordinates (Ao,{-4m})- This situation contrasts with that 
described in Eq. H4.21|l . where the total Hamiltonian is indeed Hr- This observation is one of the main motivations for 
the construction of the Fock representation (JF(7Y), {6„, fej^}) that we have presented, representation where the total 
dynamics provided by Eq. (|4.21ll is implemented in a natural way. 

(2) Let us analyze now the regime of times T € [T, oo), with T large enough so that d(\n\T) can be neglected with 
respect to the unity in an asymptotic approximation. We will then have d(|n|T) « and c(|n|T) « exp[i(5|„| (T)] 
where, for each nonzero integer the phase 5\n\{T) is some smooth real function of T. Therefore, for the nonzero 

modes, the solutions Gn {t, 9) considered in Sec. V behave at leading order as 



Gr(M)-/|e»^e-^M(^)77o(|n|t) 



(7.4) 
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for T > T. From Hankel's asymptotic expansion (|6.1|) of Hq, we also get that, for large values of t, 

G^P{t,9) « ,y|„|(T)e-^(l"l*-"^), (7.5) 
r7|„|(T) ^^e-*i"i(^)e-/4. (7.6) 

Hence, in the asymptotic region of large times, solutions (|7.5|l approach the (nonzero) "orthonormal" mode solutions 
for a free massless scalar field propagating in the static background ~ R"*" x T^, g^'^^) [see Eq. H3.2() ]. From Eqs. 
(|5.5(l and (|7.5|l it follows that, in the limit in which the system becomes massless, Jt approaches the (counterpart of 
the) Poincare-invariant complex structure of Minkowski spacetime, namely Jm = —i—£t£t)~^^'^£t, where £t is the 
Lie derivative along t" := (d/dt)'^. In fact, this is not unexpected: for asymptotically large values oft, the Hamiltonian 
(|4.9|) describes (in the sector of nonzero modes) a collection of harmonic oscillators with frequency w„ = as we 
have already seen. 

(3) In Sees. V and VI, we have formulated the quantum time evolution in the Heisenberg picture. For com- 
pleteness, we will now discuss the Schrodinger picture. In this picture, the evolution is attained by implement- 
ing the Bogoliubov transformation H4.21(l on the "one-particle" Hilbert space Ti [l^. Namely, the initial state 
("""(^o) = J2nez^ri{t,0)bn{to) at t = to w ill evolve to the final state C^itf) = J2n<£z'^riit,9)bn{tf) at t = tf, 
with bn(to) and bn{tf) related by Eq. 14.21|l . This transformation defines a pair of bounded linear maps a : 7i — s- 7i 
and (3 : TC ^ H (recall that is the complex conjugate space), 

oo 

a-C+(io) - Gn{t,e)an{t,to)bn{to), (7.7) 

n— — oo 
oo 

/3-C+(io) = E G*_^{t,e)f3:{t,to)bM- (7.8) 

n— — OO 

Unitary implementability is possible if and only if the operator /5 is Hilbert-Schmidt, that is, if and only if 

oo oo 

tr(/3t/3) = J2 {Gn,f}^P-Gn)H= J2 {P-Gn,f3-Gn)n 

n— — OO n— — oo 

oo oo 

= {p:g*_^,p:g*_j^= y \Pn\^<^- (7-9) 

n— — 00 n— — 00 

As we have shown, this is in fact the case. 

It should be stressed that the maps a and (3 define the unitary map U that implements the dynamics on T(H). 
Namely, considering the standard annihilation operator associated to C''^(*o) i-e. the smeared annihilation operator 



b (C+(to)) , we know that U is defined -up to a phase- by [22, U. 



U{tfMHC+{h))U\tfM) = &(a • C+(io)) - b^(/? • C+(io)). (7.10) 

Since the antilinear part (3 of the Bogoliubov transformation is not null, the vacuum [the state |0,<o) G ^0^) 
annihilated by b(?7) for all 77 e 7i] does not remain invariant under the action oiU. That is, Z^(t/, io)|0, to) will not 
be annihilated by b{fj) for all -q ^ Ti. Note that Eq. (|7.1U|I is just the smeared version of Ea. H5.8|) . In terms of the 
considered UF of Fock representations, the operator (|7.10() can be viewed as the annihilation operator in the tj Fock 
space (that associated with the Cauchy surface t/). Of course, the vacuum state of the tj Fock representation will 
not coincide with |0,to)i but rather be given by |0, t/) = Z/^(t/, to)|0, to)- 

Some additional comments may be worth at this stage. As we have said, the state Z^(t/, to)|0, to) G ^{'H) can be 
identified as the vacuum of the t/ Fock space. Thus, the evolution map W(t/,to) in T{H) can also be viewed as the 
unitary map relating the to and t/ Fock representations. On the other hand, in order to determine the evolution in 
the dense subspace of states with a finite number of particles, one only needs to know how the "n— particle" states 
evolve, and this in turn becomes completely fixed by specifying how the vacuum and the creation operators change 
in time: given a "n— particle" state |n^) = i'^(Ci^)b^(C2^) • . . b^'(C^)|0, to) [and abbreviating Zi(t/, to) to U], one has 

u\n^) ^m\QtV^m\c,+)u^u . . .m\C)^^u\^M)- (7.ii) 

Therefore, the adjoint of Eq. (|7.1Q|I and the corresponding relation between vacua that provides W(t/, to)|0, to) 
determines indeed the evolution in JF(7i). 
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Similar comments apply to the to Fock representation constructed without smearing operators, i.e., the Fock 
representation {J^{TL), {6„(to), &|i(io)})- One only has to replace the smeared operators b with the 6„'s and notice that 
the unitary operator is now defined by Ea. H5.8|l . Let us explain this point in more detail, for the sake of clarity. In 
the Fock representation (jF(7i), {6„(to), &l(io)})i the evolution operator is defined by 

C/(t/,io)&n(io)f/^(i/,to) =a«(i/,to)6«(io)+/3„(i/,to)&l„(io)- (7.12) 

The corresponding vacuum state |0,fo) (characterized by the conditions 6„(io)|0,to) ~ for all n) evolves 
to the state |0,t/) = C/(t/,io)|0,to), which corresponds in turn to the vacuum of the tf Fock representation 
(J-'('Wtj.), {6„(t/), 6jj(t/)}), where the annihilation operator bn{tf) is given by Eq. H7.12() . Since /3„(t/,to) does not 
vanish for tf ^ to, |0, to) and |0, tf) do not simply differ by a phase. The explicit relation between these vacua will be 
presented below. 

(4) Let us analyze now the issue of "particle" production, focusing our discussion on the nonzero modes. Among the 
UF of Fock representations ^{T{T-It), {bn{T),bl^{T)})^ ^, let us consider the T = to and T = tf representations. 

The expectation value of the number operator at time tf, namely N{tf) = '^n^Qbl^{tf)bn{tf), in the vacuum state 
at time to, \0,to), is given by 

oo oo 

{OMN{tf)\0,to) = \Mtf,to)f^2Y,\f3^{tf,to)f. (7.13) 

n— — oo,n^0 m— 1 

This expectation value is different from zero, but also bounded from above, because the sequence {f3m{tf,to)} is square 
summable for all times to-,tf > 0. As we have seen, for asymptotically large values of and tf we can neglect the 
value ofd{\n\to) and d{\n\tf), so that /3„(i/, to) ~ and an{tf,to) ~ exp{i[S\n\{tf) — 5\n\{to)])- Then, in the asymptotic 
region, &„(t/) and 6„(to) differ only by the phase exp{i[6\n\{t f) — (5|„|(to)]) [i-e., (Jt/ — J) ~ 0], and N{to) ~ N{tf). 
That is, the "particle" production decreases as to and tf grow and, consequently, |0,t/) « |0,to)- 

Actually, the evolution of the vacuum can be straightforwardly calculated by remembering that the vacuum at time 
T is characterized (up to a phase) as the unit state annihilated by all of the operators &„(T). From the evolution of 
these operators, it is then not difficult to see that the relation between the studied vacua is 



|0,i/) = Fexp 



Y,^mitf,to)bl{to)blM 



|0,to), (7.14) 



where \m{tf,to) ■= l3m{tf,to)lQ-m{tf,to) is the ratio between Bogoliubov coefficients and F is a normalization factor. 
Demanding that the vacua have unit norm, we obtain 



1^1= n Vl-|A™(</,io)P- (7.15) 

m— 1 

To the best of our knowledge expressions (of the form of) 1)7.1411 and H7.15|l were first derived, in a cosmological 
context, almost forty years ago [23. Notice that |Fp = |(0, to|0, t/) p. That is, |Fp is the probability of finding no 
"particles" [corresponding to bm{tf)] at time tf, provided that the state contains no "particles" at time to. Therefore, 
|Am(i/,io)P gives the probability of observing a nonzero number of "particles" in the mode m (or — m) at time tf. 
For asymptotically large values of to and t/, the latter probability tends to zero (so that \F\'^ approaches the unity) 
and the vacua |0,to) and |0, t/) become indistinguishable. Hence, the vacuum tends asymptotically to be stable. 

Let us remark that, as an immediate consequence of Eq. H7.14|l . "particles" are created in pairs. Besides, it should be 
emphasized that "particle" production is due to the fact that N{t) does not commute with the reduced Hamiltonian, 

00 

": Rr N{t)\ = /'("M) (bm{t)b-ra{t) - blitfbl^it)) . (7.16) 

m— 1 

However, since P{m,t) = l/(8mt^), the commutator vanishes in the asymptotic limit of large times when the system 
becomes massless, i.e. when the time dependent potential equals zero. So, the theory becomes "free" and the quantum 
representation approaches (the counterpart of) the Poincare-invariant one. 
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VIII. CONCLUSION AND FURTHER COMMENTS 

Let us summarize our results, discuss some consequences of the introduced quantization and compare it with 
previous ones. The first observation that is worth emphasizing is that, even if the field (,{to;d) evolves unitarily to 
£_{t; 9), the (explicitly time dependent) formal operator (f>{t; 9) := ^(t; 9)/\/t, which was regarded as the basic field for 
the quantum model in Ref. Q, does not display a unitary evolution. In other words, 0(io; ^^nd 4>{t] ^) are not related 
by means of a unitary operator. The choice of fundamental field plays therefore a decisive role in the construction 
of a satisfactory quantization (together with the subsequent choice of annihilation and creation operators, i.e. the 
complex structure J). To arrive to the new field-parametrization, we have benefitted from the freedom available to 
introduce a time dependent canonical transformation and redistribute in that way the time dependence in an implicit 
part, whose evolution is generated by the corresponding reduced Hamiltonian of the model, and an explicit part (the 
factor for (f) in our case), whose variation does not necessarily have to be described by a unitary transformation. 
We notice that, in systems like the symmetry reduced Gowdy model where the Hamiltonian depends explicitly on 
time, it is natural to take into account the possibility of performing canonical transformations that vary with time. It 
is also worth pointing out that this system exhibits what seems to be a general feature of quantum field systems, in 
the sense that generic linear canonical transformations do not become unitarily implemented in the quantum theory 
(see for instance and p^l. 

On the other hand, we have seen that the vacuum of the quantum theory proposed for ^ is not left invariant by the 
time evolution. As a consequence, there is some "particle" production by the vacuum, which certainly attenuates as 
t becomes large, but never becomes strictly zero except in the limit in which the time dependent potential vanishes, 
i.e. at infinitely large times. However, we have seen in detail that this poses no problem for unitarity and that the 
quantum theory is perfectly consistent, with a well defined evolution that is compatible with the standard probabilistic 
interpretation of quantum mechanics. 

The fact that particles can be created (in pairs) in expanding universes was realized already in the late sixties 
|27l. Since then, particle creation has been extensively discussed and studied in diverse contexts in QFT (see e.g. 
|5l l25l 113), leading to remarkable results as the so-called Hawking and Unruh [s^l effects. In general, these 
results rest on the analysis of the (specific) Bogoliubov transformation that relates the canonical operators between 
the in and out states. In the Gowdy cosmological model, as we have seen, the positive and negative frequency 
parts of the basic field ^ become mixed during evolution owing to the time dependent potential V{£) — ^ j (4i^). The 
canonical annihilation and creation operators associated with out Fock states, at time Tout, are linear combinations 
of those associated with in Fock states, at an earlier time Tin. From the unitary implementability of the dynamics, it 
follows that every in state with a finite number of "particles" evolves to an out state which also has a finite (although 
possibly different) number of "particles" . 

This result applies even when Tout tends to infinity. In that limit, a neat particle interpretation is available 
for the out Fock states. Indeed, finiteness in the number of "particles" is a simple consequence of the fact that 
limr^^t^oo |/3ra(Tout, Tin)| = |dra(Tin)| and that dnify is square summable for all t > 0. On the other hand, the 
normalized om< modes G'n°°^(t,^) :— HmT<,ut^oo exp [j(5|„|Tout — «7r/4] gI"^""'-* (i, 0) (see Sec. [VII|l behave like the positive 
frequency modes of Minkowski spacetime (except for the different background topology) in the limit where the 
system becomes free. Therefore, states in the out Fock space .Fout('Woo), which is constructed from the space of 
solutions and the natural complex structure defined by the modes ci""' , admit a natural particle interpretation in the 
asymptotic future. That is, .Fout('^oo) can be asymptotically identified with the standard Fock representation of a free 
massless scalar field propagating in a Minkowski spacetime, so that well defined asymptotic notions of vacuum and 
particles arise. Clearly, in the limit where the system is invariant under time translations, ambiguities in the particle 
interpretation are avoided. But, furthermore, an approximate adiabatic notion of particles [Slj can be introduced for 
each finite, sufficiently large value of T, because in the asymptotic regime the potential V varies then very slowly in 
time. Thus, for large Tin and Tout, a notion of particles with a conventional interpretation is available. One should, 
of course, keep in mind that the notion of "particle" we refer to is in the sense of QFT on curved space, where the 
existence of a well defined notion of particles refers to actual particles (as registered in detectors). In the present 
quantum gravity system, even when its degrees of freedom are captured by the scalar field, the "particles" associated 
with this field are far from having a clear interpretation in terms of geometrical objects or "quanta" . 

The Gowdy model was reduced in Ref. @| to a free massless scalar field on a flat, but time dependent background, 
subject to a global constraint. We know that a conventional quantization of this field does not allow one to represent 
the evolution by means of a unitary transformation. By a field redefinition, which involves the time parameter, we have 
mapped the system into a scalar field ^ satisfying a "Klein- Gordon" type equation in a background which is flat and 
time independent (like 3-dimcnsional Minkowski spacetime, apart from the topology), although in presence of a time 
dependent potential. The natural quantization of this new field that we have presented in this work provides a theory 
with a unitary dynamics. The relation between this theory and the natural time-translation invariant quantization that 
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is available in this background (namely, the analog of the Poincare invariant quantization in the considered spacetime) 
becomes manifest in the asymptotic regime where the time-translation symmetry is recovered and a preferred vacuum 
with that symmetry can be selected. More precisely, for asymptotically large values of T, we have shown that our 
complex structure Jt approaches the "Poincare-invariant" one in the limit in which the system becomes massless 
and hence invariant under time-translations. Although this result suggests the appealing possibility that there exists 
a connection between unitary implementability and asymptotic symmetries, further research is needed to elucidate 
what kind of physical requirements lead to acceptable quantizations (an equivalent of the Hadamard condition in 
QFT in curved space pill. 

In order to arrive at a quantum description of the linearly polarized Gowdy model one needs inputs at both the 
classical and quantum levels. Classically, there are two important inputs. One is the choice of deparametrization, 
namely the choice of a fictitious time. The other is the field-parametrization adopted for the metric. We have seen 
that, once a choice of time gauge has been selected, the freedom in the choice of basic field can be understood as 
that in performing time dependent canonical transformations in the system. Besides, the order in which the two 
previously mentioned choices are made in our case is irrelevant, because it does not affect the final outcome, as we 
show in the Appendix. In the quantum part, on the other hand, the cornerstone is the choice of a complex structure 
(still a classical construct on phase space), that determines the vacuum and the structure of the quantum theory. 
In this respect, our choice of complex structure J was in some sense natural, guided both by previous proposals for 
the quantization and by the symmetries of the system in the asymptotic region of large times. It would be very 
interesting to determine all other possible quantum representations of the canonical commutation relations for our 
specific field-parametrization (and time gauge) that permit a unitary implementation of the dynamics, together with 
some additional physical requirements, and elucidate whether they are all unitarily equivalent. If this were so, the 
quantum theory presented here would be essentially unique, once the choice of internal time and fundamental field 
has been fixed. This issue will be the subject of a future investigation |3^ . 

Finally, let us present some general comments on the validity of our results in the more general context of quantum 
gravity. Quantum gravity, both in its full glory and in reduced (midi-superspace) models, suffers from the celebrated 
problem of time. Roughly speaking, this means that there is no fundamental notion of time (even classically) and 
that this has implications for the usual probabilistic interpretation in the quantum theory. The Gowdy models are 
not free from such a problem. There are basically two different approaches, resulting from the two different ways 
of quantizing constrained systems: quantize first and then reduce (Dirac) or reduce first an then quantize (reduced 
phase space). The procedure that has been followed here and in Refs. JXlj JJJ for the Gowdy model, although 
closer to the second option, is a mixture of both approaches: one reduces the system classically via gauge fixing and 
deparametrization, but keeps a global constraint, which is dealt with in the quantum theory. For the problem of time, 
one is choosing an internal time i, via the deparametrization procedure, that takes the role of an "external parameter" 
in ordinary quantum theory. Of course, the parameter t is not the physical time even classically, but provides us with 
the familiar framework of quantum theory to answer real physical questions. 

From the viewpoint of canonical quantum theory, where the theory is defined over an abstract 3-manifold S, the 
natural picture for the quantum description of gravity is the Heisenberg picture. A quantum state \^) of the system 
is defined on S (that should not be thought of as a "constant-time slice" since there is no time and no spacetime to 
embed it), and observables are of the Heisenberg type, namely evolving constants of motion p^.ll6j. To be precise, we 
have for the Gowdy model a family of operators {Oi(t)}, one observable for each value of t. One can, of course, define 
these operators and relate them by means of the "evolution operator" U{t, to) which is, as we have shown, perfectly 
well defined. The Heisenberg picture has to be contrasted with the Schrodinger one that in full quantum gravity (in 
the Dirac approach) is simply not defined, since there is no notion of embedding of the hypersurface S on a spacetime, 
much less the notion of "time evolution" . In our case, however, since we have defined a notion of (internal) time, 
we are free to try and construct both the Schrodinger and the Heisenberg pictures. As we have shown, these two 
pictures are well defined in our model. Finally, let us end this note by pointing out that in order to make full justice 
to the quantum geometry description given by our choice of quantization, one would need to analyze the behavior 
of (quantum) metric objects that provide a description of the quantum geometry as in Ref. j^, and the observables 
recently introduced in Ref. |33]. We shall leave that analysis for future research. 
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APPENDIX A: GAUGE FIXING IN THE NEW FIELD-PARAMETRIZATION 



In this appendix, we explicitly show that the canonical transformation H3.9|l amounts to a field-reparametrization 
of the metric of the Gowdy model which commutes with the process of gauge fixing p^ . In particular, the gauge 
choice is not modified. 

Let us start with the 3 + 1 decomposition of the metric for the polarized Gowdy spacetimes after fixing the 
gauge corresponding to diffeomorphisms in the direction of the coordinates a E and 5 E , with the two axial 
Killing vector fields identified with d^r and ds |12] '■ 

ds^ = -N^dt^ + hgg[de + N'^dtf + h„„da^ + hssdS"^. (A.l) 

Here, N is the lapse function, is the 6'-component of the shift vector, and hij is the induced spatial metric. To 
arrive at this expression, we have employed the fact that the two Killing vector fields are hypersurface orthogonal, so 
that haS must vanish. Besides, the presence of the Killing symmetries implies that all metric functions depend only 
on 6 £ and on the time coordinate t, which we choose to be positive. 

Instead of adopting the same field-parametrization as in Ref. for the induced metric, namely 

hee = e'^"'^, h^a = e~'^T'^, hss = e'^, (A. 2) 

we now introduce an alternative parametrization in terms of a new set of fields {Q"} ■= {^,7,''"}, 

/,,, = e^-(«/v^)-«^/(4-), /^.. = e-^/v^r^ hss^e^/^. (A.3) 
With this new field-parametrization, we obtain 

ds^ = e^-(«/v/^)-«V(4r) (^r^N^dt^ + [de + N'dtf) + e'^'^ (r^do^ + e^^l^db"^ , (A.4) 

where N N j \fh is the densitized lapse function and h is the determinant of the induced metric. 

Regarding the change from the components of hij to the set {Q"} as a point transformation, it is straightforward 
to find momenta Pa canonically conjugate to in terms of those for the induced metric [T^ . In this way, one arrives 
at the following Einstein-Hilbert action in Hamiltonian form: 



/."*/ 



dB 



PrT + PjJ + P^i - {NC + N^Cg 



(A.5) 



where the momentum and (densitized) Hamiltonian constraints adopt the respective expressions 

Cg = P^r' + B^f + F^e'-SP;, (A.6) 
... , r 



C = -Pi + i^P^-TPrP^i 



2 



2 « 8t ^ ' ' 2 
On the other hand, a comparison between the field-parametrizations ljA.2|) and ljA.3|l shows 



(A.7) 



7 = 7-|^, V=^, (A.8) 

with the same field r in both cases. This point transformation leads then to the following relations between the 
corresponding canonical momenta: 

where we have called Pr the momentum conjugate to r in the old parametrization, to distinguish it from the new 
one, Pr- 

In order to deparametrize the model and fix (almost all of) the remaining gauge freedom, we must impose additional 
conditions that, together with the constraints ljA.6|) and ljA.7|l . form a set of second class constraints allowing the 
reduction of the system. In Ref. (l2l| . the conditions imposed were gi := P^ + p ~ and g2 ■= t — tp = 0, where 
p = — § d9Pj/{2TT) is a constant of motion for the model. Note that the only canonical variables that appear in these 



19 



conditions are and r. With our change of metric fields, r is not modified and Pj becomes P^. Therefore, the 
gauge fixing selected to arrive at the usual description of the Gowdy model, and that we choose to impose also in the 
new field-parametrization, is 



gi := Pj + p ^ 0, g2 = T - tp = 0. 



(A.IO) 



An analysis along the lines discussed in Ref. ^3 shows that the gauge fixing is well posed provided that p ^ 0. As 
explained in the main text, we restrict all considerations to the sector of positive p. The compatibility of the gauge 
fixing with the dynamics sets N = \/{pt), whereas can be any fimction of t. Although the shift vector is not 
entirely determined, the allowed functional form is such that its contribution to the metric can be absorbed by means 
of a redefinition of the coordinate 9. 

The momentum constraint Cg = Q together with the gauge fixing conditions imply that 



Pi' = P^^'- 



(A.ll) 



Since p > 0, this relation determines the function 7, except for its zero mode. Given the periodicity of the system in 
6*, Eq. (|A.lip also supplies the homogenous constraint that remains on the system. 



Co 



0. 



As a result of the commented gauge fixing, one finally obtains a reduced system with spacetime metric 



Q . ^ 1 

7 — 1 y 

27rn ^-^ 2Tmp 



d9^ 



e-^/^fp^da^ + e^^^dS^ 



d0e™(«-«)Pf£' + — ([d9 



Anp J 



Pi + (O' + 



(A.12) 

(A.13) 
(A. 14) 



where Q is the configuration variable canonical to P = Inp G M. This metric coincides in fact with that obtained 
from Eq. (|2.1|) (expressed in terms of the field (p) when the time dependent canonical transformation (|3.9|) is applied 
directly in the reduced model. 

In addition, the action that one obtains after the gauge fixing procedure is (modulo a spurious boundary term) 



dtlPQ 



d9 



/\2 



1 

4^2 



(A.15) 



So, the reduced Hamiltonian is precisely that deduced in Eq. H3.12|l . Therefore, as we wanted to show, the field- 
reparametrization commutes with the gauge fixing and reduction procedure. 
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